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Abstract — An accurate channel estimation is crucial for the 
novel time domain synchronous orthogonal frequency-division 
multiplexing (TDS-OFDM) scheme in which pseudo noise (PN) 
sequences serve as both guard intervals (GI) for OFDM data 
symbols and training sequences for synchronization/channel 
estimation. This paper studies the channel estimation method 
based on the cross-correlation of PN sequences. A theoretical 
analysis of this estimator is conducted and several improved 
estimators are then proposed to reduce the estimation error floor 
encountered by the PN-correlation-based estimator. It is shown 
through mathematical derivations and simulations that the new 
estimators approach or even achieve the Cramer-Rao bound. 

I. Introduction 

In the recently proposed time domain synchronous orthogo- 
nal frequency-division multiplexing (TDS-OFDM) scheme jl], 
the classical cyclic prefix (CP) conventionally used in OFDM 
is replaced by a known pseudo noise sequence (simply termed 
as PN hereafter) which is reused as training sequence for 
channel estimation and synchronization. Consequently, TDS- 
OFDM combines the guard interval (GI) and the training 
symbols and does not need any additional pilots in the fre- 
quency domain, thereby achieving a higher spectral efficiency 
than CP-OFDM. TDS-OFDM has been adopted by the novel 
Chinese digital television broadcasting standard-DTMB IS). 

In TDS-OFDM, a channel estimate is needed to separate the 
PN from the OFDM data part at the receiver Its accuracy is 
crucial for the demodulation process to avoid any residual PN 
components in the received signal. Hence, channel estimation 
plays a prominent part in TDS-OFDM performance and needs 
to be carefully studied. [3] investigates several channel estima- 
tion techniques based on known sequences but did not exploit 
the property of the PN. JT] uses the cross-correlation results 
of the transmitted and received PNs as the channel estimates 
ignoring the fact that the correlation of the PN is not a 
perfect impulse function, which introduces mutual interference 
between different channel paths in the estimation result. To 
overcome this problem, (|4] suggests to remove interference 
components of several strongest paths and fSl proposes to 
iteratively detect the significant paths in the correlation results. 
Both of them did not thoroughly eliminate the interference 
lying in the PN correlation. Finally, [6] proposes a least square 
(LS) channel estimator based on the property of maximum 
length sequence (m-sequence). This solution however gives a 
suboptimal performance as explained in the sequel. 
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Fig. 1. Structure of the GI specified in the DTMB standaid. 



In this paper, we investigate the performance of the PN- 
correlation-based channel estimator, and, more importantly, 
propose several improved estimators that reduce or even elim- 
inate the estimation error floor resulting from the interference 
term issued from the PN correlation function. In the following, 
section II introduces the PN based channel estimation and 
gives a theoretical analysis of its performance. In section III, 
three different improved estimators are proposed and their 
respective mean square errors (MSE) are derived. Finally, their 
performance is compared in section IV through simulations 
over different channel conditions. 

II. PN-CORRELATION-BASED CHANNEL ESTIMATOR 

A. System Model 

In the DTMB system, the i^-length GI consists of an A^- 
length m-sequence as well as its pre- and post- circular 
extensions [21 . Since any circular shift of an m-sequence is 
itself an m-sequence ||7|, the GI can also be treated as another 
A^-length m-sequence, denoted by p = [pi,p2, • ■ • iVnY with 
(•)^ standing for the matrix transpose, and an A'pp -length CP. 
That is, if this CP is longer than the channel delay spread, 
the CP absorbs the channel time dispersions and the A^-length 
m-sequence is ISI-free. The received PN can be written as: 

d =Hp + w = Ph + w, (1) 

where d = \d\,di, . . . ,df^Y ^ ^ = \w\,wi, . . . .w^i^ and 
h = [ft-i, h,2i . . . , ft-L, 0, . . . , 0]^ are the received PN, additive 
white Gaussian noise (AWGN), and channel impulse response 
(CIR), respectively. The channel is modeled as an L*-order 
finite impulse response (FIR) filter with = E[|/i;p] the 



average power of the channel tap. The power of the 
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channel paths is normalized such that J2i=i^ 
and P are the N x N circulant matrices with first rows 
[/ii,0,...,0,/iL,...,/i2] and [pi,pN,PN-i, ■ ■ ■ ,P2], respec- 
tively. The second equality in ([T]i uses the commutativity of 
the convolution. 



B. PN Correlation-based Estimator 

Channel estimation can be simply obtained by performing 
time domain correlation of transmitted and received PN se- 
quences. Recall the circular autocorrelation property of the 
m- sequence Q: 
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where (•)* is the conjugate of complex number and [ Jjv 
denotes modulo-A^ operation. When the PN is sufficiently 
long, i.e., — w 0, the autocorrelation function approaches 
the Kronecker delta function. The CIR can be extracted from 
the received PN using this delta function. The estimator is: 
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where C = is an TV x circulant matrix with first row 
\pitP2, ■ ■ ■ iPn\ and represents the circular correlation of the 
PN. {■)^ is the Hermitian transpose. Using the correlation 
property ([2), it yields: 
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The estimation error is = h — h = (Q — I^vjh 
where Ijv is an TV x TV identity matrix. Assuming that the 
channel taps are uncorrected, the MSE of the estimate is: 

= 1 Tr (E [ (Q - lAr)hh«(Q - I^f] )+^Tr (e [Cww« C^ 
=lTr ((Q - I^)A(Q - + ^Tr (cC«) , (5) 



interference noise 

where Tr(-) is the trace of a matrix, A is a diagonal matrix 
of size TV whose first L elements of the main diagonal is 
[al,a2, . . . , , and the rest elements are all O's. The estima- 
tion error is composed of two parts: the interference resulting 
from the correlation of PN sequences and the noise. It can 
be found that the interference comes from the contributions 
of the off-diagonal elements in Q and vanishes if Q = In- 
The estimator is asymptotically unbiased as TV — > oo. More 



concretely, the interference is computed as: 
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Therefore, the MSE of the estimator is finally: 

a| TV^ 
TV TV3 ■ 

Eventually, the first term of the MSE expression is proportional 
to the noise variance, while the second term is only determined 
by the length of the PN. In other words, it produces an 
estimation error floor with an MSE level of (TV — 1)/TV'^. 
This error floor appears when (TV — 1)/TV'^ > ct^/TV, i.e., 
SNR(dB) > lOlogio (TVV(^ - !))■ For example, the error 
floor appears when the SNR is greater than 24.1 dB and 
27.1 dB, given the 255-length and 511-length PNs specified 
in ID, respectively. 

The Cramer-Rao bound of training sequence based channel 
estimation with a length equal to the training sequence is lH): 
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Comparing O and (|8]l, it can be found that the correlation 
based estimator approaches the Cramer-Rao bound at low 
SNR (i.e. cr^ is large), but suffers an estimation error floor at 
high SNR. Therefore, we go in for some improved estimators 
aiming at reducing this error floor 

III. Improved Estimators with Reduced Error 
Floor 

A. Method 1: Multiplying by Inverse of Matrix Q 

From the analysis in the last section, the estimation error 
floor comes from the fact that Q is not a perfect identity 
matrix. Therefore, a straightforward solution is to perform a 
linear transformation Q, such that fJQ = Iat. Since Q is known 
and always full rank for a given m-sequence, O = Q^. A 
new estimator is obtained by left multiplying the correlation- 
based estimator (|3]l by Q^: 



hi =Q-ih 



(9) 



which leads to an LS estimator that is in some extend similar to 
that proposed in IS). The estimation error is = iQ^^Cw. 
The MSE of the estimator is: 
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As Q is a circulant matrix and all its elements are known for a 
given m-sequence, its inverse is easily obtained and does not 
need complex computations ||9l : 



where 



2N 



b b 



■, and b = 



N 



(11) 



N+r N+i 

Replacing (fTTI ) into ( fTOb . the MSE becomes: 
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Comparing (fT2] i and (|7]i, it can be found that the estimation 
error floor is removed by left multiplying by Q~^. In the 
meantime, the MSE is however approximately twice as much 
as the Cramer-Rao bound due to a noise power increase. 

B. Method 2: Multiplying by Inverse of Truncated Matrix Q 

Suppose the length of the CIR L is perfectly known, the 
estimate (O is truncated to L-length. The estimator becomes: 
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where h is the L-length vector of the real CIR. T is an L x iV 
matrix whose left L x L submatrix is an identity matrix and 
rest parts are all O's, which represents deleting the last N — L 
elements of an A^-length vector. Q = TQT^ is an L x L 
circulant matrix that contains the first L rows and L columns 
of matrix Q and is still full rank. Similarly to (|9]l, we can left 
multiply estimator ( fT3] l by the inverse matrix Q to obtain: 

h2 =Q"^h = h+^Q"^TCw. (14) 

This new estimator is still unbiased. The related estimation 
error is = -^Q TCw which leads to an MSE given by: 

1^ /^r. ,H- 



= §Tr(Q-Q(Q-V).||Tr(Q-).(15) 

The inverse matrix Q ^ is the L x L circulant matrix with a 
similar form as (fTTI . The elements a and b are replaced by 
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N'^ + 2N - NL- L + l' 
The MSE of the estimation is finally: 

N-L + 2 

^'^^ ^ N'^ + 2N -NL-L + l 



(16) 



Comparing (fTSl i and ( fT6] l, it can be found that the estimation 
MSE is reduced thanks to the truncation process. In the 
extreme case when N = L, the estimator ( fT4] i turns to (|9). 



When the channel length L is known, the Cramer-Rao 
bound of a training sequence based channel estimation is 
computed by replacing matrix P by P = PT"^ in ([8]): 
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m + 2N - NL- L + l 
This demonstrates that the proposed estimator (ITSt achieves 
the Cramer-Rao bound. 

C. Method 3: Subtracting Interference 

Another improved estimator can be obtained by subtracting 
the contribution of the interference from the correlation-based 
estimator. 

More precisely, suppose the CIR length L is known, we can 
rewrite the truncated CIR estimate (fTsl i as: 

1 
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where A is the Lx L matrix that contains all the off-diagonal 
elements of Q and the main diagonal elements of A are O's, 
h is the L-length vector of the real CIR. Considering that 
A is known for a given m-sequence, we propose to use the 
estimated h to reduce the interference. An estimator is thus: 
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The estimation error is: = — Ah + -i-TCw — -^ATCw. 
RecaUing that Q = -^TCC^T^, the MSE of the estimate is: 

= iTr^AA(A)«+^(Q+AQA''-2K{QA''})) (21) 

where A is an L x L diagonal matrix with diagonal elements 
[af , ■ ■ ■ , aW, 3?{-} standing for the real part of a complex 
number. The first term of the MSE expression is: 

AA(A)^ = 
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The third term of the MSE expression is: 



^AQA^ 



Ar4 



The fourth term of the MSE expression is: 
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Fig. 2. MSB performance of different estimators with v = 420, N = 255, 
7V„„ = 165 in tlie TU-6 channel. 



TABLE I 

Complexities of the PN-based channel estimation methods. 



Estimator 


Complexity 


Circular con'elation-based estimator (|3] 


0(Af2) or 
O(Ar-logAf) 


Multiplying matrix inverse Af-length ^9) 
Multiplying matrix inverse L-length )14t 
Subtracting inteiference ||19I 


O(N^) 



Eventually, putting (|22li, (|23ll and (l24]i into (EB yields: 

_ N^ + {L-l){2-L-N) 2 {L-1){L^ -3L + 3) 

(25) 

Comparing (|25]l with (Q, in low SNR region, the MSE 
approaches the Cramer-Rao bound, while, in high SNR region, 
the estimation error floor is reduced by a ratio of approxi- 
mately (^)'. 

IV. Complexity Analysis 

In order to compare the computational complexity of each 
estimation method, we investigate the required basic oper- 
ations, i.e. multiplication, additions and FFT's etc. Table H] 
shows the complexity of the basic correlation-based estimator 
(O and the additional complexities required by the improved 
methods (|9|, (O and ( fT9] l. 

The A^-length circular convolution takes N multiplications 
and TV — 1 additions for each delay. That is N'^ multiplications 
and X (A^ — 1) additions for an A^-length CIR estimate. 
Therefore, the complexity of the circular correlation-based es- 
timator ^ is 0{N^). Considering that the circular convolution 
can also be computed by using the FFT, the computational 
complexity can be reduced to 0{N ■ log A^). 

As far as the improved estimators are concerned, additional 
complexities are needed by the estimation refinement process. 
For instance, in the estimator (|9]l, the matrix Q^^ depends only 
on the given PN sequence and can thus be pre-computed and 



Fig. 3. MSE performance of different estimators with u = 945, N = 511, 
N^„ = 434 in the TU-6 channel. 



TABLE II 

Profile of the TU-6 and HT channels 



Channel 


Tapl 


Tap2 


Tap3 


Tap4 


Taps 


Tap6 


TU 


Delay (/is) 





0.2 


0.5 


1.6 


2.3 


5.0 


Power (dB) 


-3 





-5 


-6 


-8 


-10 


HT 


Delay (/is) 





0.2 


0.4 


0.6 


15.0 


17.2 


Power (dB) 





-2 


-4 


-7 


-6 


-12 



stored 0. The matrix multiplication needs A^^ multiplications 
and A^ x (A^ — 1) additions. Therefore, the additional com- 
plexity is 0{N^). In the meantime, the matrix multiplication 
can be carried out in a reduced length L as done in estimator 
(fT4l) . For the same reason as in the pervious situation, the 
computation of Q does not need any additional effort. The 
corresponding additional complexity is therefore reduced to 
0{L^). For the estimator ([19}, the refinement of each channel 
tap needs {L — 1) additions. That is, L x (L — l) additions are 
required for the whole CIR estimate, which corresponds to a 
computational complexity of 0{L^). 

V. Simulation Results 

The simulation parameters are selected according to the 
specifications of the DTMB standard 121 where the sampling 
period is 1/7.56 /xs. The PN sequences are generated using 
the maximal linear feedback shift registers specified by the 
standard as well. The Typical Urban with six paths (TU-6) and 
Hilly Terrain (HT) channels specified in [TTol are used in the 
evaluation. The power delay profiles of the two channels are 
given in Table [III The maximum delays of the TU-6 channel 
and HT channel are 5 /is and 17.2 /is which correspond to i = 
38 and L = 130 samples of the DTMB system, respectively. 
The performance of the different estimators is investigated 

'in fact, the matrix Q^' is quite structured. There are only two values - 
diagonal and off-diagonal elements. Hence, it is not necessary to store all the 
elements of the matrix. Instead, it is smarter to record the two values only. 
The cost of the storage is therefore negligible. 
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Fig. 4. MSB performance of different estimators with v = 420, N - 
Nf,„ = 165 in tlie HT channel. 
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Fig. 5. MSE performance of different estimators with u - 
AT = 434 in the HT channel. 



945, N = 511, 



with different PN lengths and channels as shown from Figure |2] 
to Figure |5] Results are obtained from 1000 realizations of 
each channel. The correlation-based estimator corresponds to 
the classical approach found in the literature [l] while others 
are the improved versions introduced in section |III]A, B and 
C. From the figures, we can have the following observations: 

1 . Methods 'improved method 1 ' and 'improved mettiod 2' com- 
pletely eliminate the estimation error floor, while method 
'improved metliod 3' reduces the error floor when the channel 
delay spread is shorter than the PN length. Multiplying the 
basic PN-correlation-based estimator by the inverse of the 
correlation matrix can create a perfect identity matrix which 
leads to estimators free of estimation error floor This can be 
seen from the performance of 'improved method 1' and 'improved 
method 2'. On the other hand, 'improved method 3' subtracts 
the interference components using the CIR estimates. It ap- 
proaches the Cramer-Rao bound in low SNR region, while the 
estimation error floor is reduced by a ratio of approximately 

in high SNR region. It indicates that the reduction is 
more notable when the channel delay is significantly small 
compared to the PN length. For instance, comparing Fig. |2] 
and Fig. IH 'improved method 3' obtains more improvement in 
the TU-6 channel than in the HT channel given the same PN. 
In addition, 'improved method 3' can achieve a lower estimation 
error floor when a longer PN is used (comparing Fig. Eland 
Fig. [51 or Fig. IH and Fig. [Sj. 

2. Method 'improved mettiod V boosts the noise variance, 
while methods 'improved metliod 2' and 'improved metliod 3' do 
not. From this point, the latter two methods are the preferred 
ones because they outperform the classical correlation based 
approach whatever the SNR. In contrast, 'improved method 1' 
presents a performance back-off due to the noise component 
power boost, which leads to worse MSE results than the 
correlation based approach at low SNR. Last but not least, 
the estimator 'improved method 2' does not have estimation error 
floor and achieves the Cramer-Rao bound. 



VI. Conclusion 

In this paper, we have investigated the PN-correlation- 
based channel estimator for TDS-OFDM. Aiming at reducing 
the estimation error floor encountered by the classical PN- 
correlation-based estimator, we have proposed three improved 
estimators which exploit the correlation property of the m- 
sequence and the knowledge of the channel length. It has been 
shown through mathematical derivations and simulations that 
the new proposed estimators approaches or even achieve the 
Cramer-Rao bound. 
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Abstract — An accurate channel estimation is crucial for the 
novel time domain synchronous orthogonal frequency-division 
multiplexing (TDS-OFDM) scheme in which pseudo noise (PN) 
sequences serve as both guard intervals (GI) for OFDM data 
symbols and training sequences for synchronization/channel 
estimation. This paper studies the channel estimation method 
based on the cross-correlation of PN sequences. A theoretical 
analysis of this estimator is conducted and several improved 
estimators are then proposed to reduce the estimation error floor 
encountered by the PN-correlation-based estimator. It is shown 
through mathematical derivations and simulations that the new 
estimators approach or even achieve the Cramer-Rao bound. 

I. Introduction 

In the recently proposed time domain synchronous orthogo- 
nal frequency-division multiplexing (TDS-OFDM) scheme [?], 
the classical cyclic prefix (CP) conventionally used in OFDM 
is replaced by a known pseudo noise sequence (simply termed 
as PN hereafter) which is reused as training sequence for 
channel estimation and synchronization. Consequently, TDS- 
OFDM combines the guard interval (GI) and the training 
symbols and does not need any additional pilots in the fre- 
quency domain, thereby achieving a higher spectral efficiency 
than CP-OFDM. TDS-OFDM has been adopted by the novel 
Chinese digital television broadcasting standard-DTMB [?]. 

In TDS-OFDM, a channel estimate is needed to separate the 
PN from the OFDM data part at the receiver Its accuracy is 
crucial for the demodulation process to avoid any residual PN 
components in the received signal. Hence, channel estimation 
plays a prominent part in TDS-OFDM performance and needs 
to be carefully studied. [?] investigates several channel estima- 
tion techniques based on known sequences but did not exploit 
the property of the PN. [?] uses the cross-correlation results 
of the transmitted and received PNs as the channel estimates 
ignoring the fact that the correlation of the PN is not a 
perfect impulse function, which introduces mutual interference 
between different channel paths in the estimation result. To 
overcome this problem, [?] suggests to remove interference 
components of several strongest paths and [?] proposes to 
iteratively detect the significant paths in the correlation results. 
Both of them did not thoroughly eliminate the interference 
lying in the PN correlation. Finally, [?] proposes a least square 
(LS) channel estimator based on the property of maximum 
length sequence (m-sequence). This solution however gives a 
suboptimal performance as explained in the sequel. 
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Fig. 1. Structure of the GI specified in the DTMB standaid. 



In this paper, we investigate the performance of the PN- 
correlation-based channel estimator, and, more importantly, 
propose several improved estimators that reduce or even elim- 
inate the estimation error floor resulting from the interference 
term issued from the PN correlation function. In the following, 
section II introduces the PN based channel estimation and 
gives a theoretical analysis of its performance. In section III, 
three different improved estimators are proposed and their 
respective mean square errors (MSE) are derived. Finally, their 
performance is compared in section IV through simulations 
over different channel conditions. 

II. PN-CORRELATION-BASED CHANNEL ESTIMATOR 

A. System Model 

In the DTMB system, the i^-length GI consists of an A^- 
length m-sequence as well as its pre- and post- circular 
extensions [?]. Since any circular shift of an m-sequence is 
itself an m-sequence [?], the GI can also be treated as another 
A^-length m-sequence, denoted by p = [pi,P2, • ■ • iVnY with 
(•)^ standing for the matrix transpose, and an A'pp -length CP. 
That is, if this CP is longer than the channel delay spread, 
the CP absorbs the channel time dispersions and the A^-length 
m-sequence is ISI-free. The received PN can be written as: 

d =Hp + w = Ph + w, (1) 

where d = \d\,di, . . . ,df^Y ^ ^ = \w\,wi, . . . .w^i^ and 
h = [ft-i, h,2i . . . , ft-L, 0, . . . , 0]^ are the received PN, additive 
white Gaussian noise (AWGN), and channel impulse response 
(CIR), respectively. The channel is modeled as an L*-order 
finite impulse response (FIR) filter with = E[|/i;p] the 



average power of the channel tap. The power of the 
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channel paths is normalized such that J2i=i^ 
and P are the N x N circulant matrices with first rows 
[/ii,0,...,0,/iL,...,/i2] and [pi,pN,PN-i, ■ ■ ■ ,P2], respec- 
tively. The second equality in ([T]i uses the commutativity of 
the convolution. 



B. PN Correlation-based Estimator 

Channel estimation can be simply obtained by performing 
time domain correlation of transmitted and received PN se- 
quences. Recall the circular autocorrelation property of the 
m-sequence [?]: 
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1 n = 

0<n< N 
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where (•)* is the conjugate of complex number and [ Jjv 
denotes modulo-A^ operation. When the PN is sufficiently 
long, i.e., — w 0, the autocorrelation function approaches 
the Kronecker delta function. The CIR can be extracted from 
the received PN using this delta function. The estimator is: 
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(3) 



where C = is an TV x circulant matrix with first row 
\pitP2, ■ ■ ■ iPn\ and represents the circular correlation of the 
PN. {■)^ is the Hermitian transpose. Using the correlation 
property ([2), it yields: 
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The estimation error is = h — h = (Q — I^vjh 
where Ijv is an TV x TV identity matrix. Assuming that the 
channel taps are uncorrected, the MSE of the estimate is: 

= 1 Tr (E [ (Q - lAr)hh«(Q - I^f] )+^Tr (e [Cww« C^ 
=lTr ((Q - I^)A(Q - + ^Tr (cC«) , (5) 



interference noise 

where Tr(-) is the trace of a matrix, A is a diagonal matrix 
of size TV whose first L elements of the main diagonal is 
[al,a2, . . . , , and the rest elements are all O's. The estima- 
tion error is composed of two parts: the interference resulting 
from the correlation of PN sequences and the noise. It can 
be found that the interference comes from the contributions 
of the off-diagonal elements in Q and vanishes if Q = In- 
The estimator is asymptotically unbiased as TV — > oo. More 



concretely, the interference is computed as: 

(Q-lAr)A(Q-lN)^ = 
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Therefore, the MSE of the estimator is finally: 

a| TV^ 
TV TV3 ■ 

Eventually, the first term of the MSE expression is proportional 
to the noise variance, while the second term is only determined 
by the length of the PN. In other words, it produces an 
estimation error floor with an MSE level of (TV — 1)/TV'^. 
This error floor appears when (TV — 1)/TV'^ > ct^/TV, i.e., 
SNR(dB) > lOlogio (TVV(^ - !))■ For example, the error 
floor appears when the SNR is greater than 24.1 dB and 
27.1 dB, given the 255-length and 511-length PNs specified 
in [?], respectively. 

The Cramer-Rao bound of training sequence based channel 
estimation with a length equal to the training sequence is [?]: 

2 2 

MSE > ^Tr ((P^P)-i) = (8) 

Comparing O and (|8]l, it can be found that the correlation 
based estimator approaches the Cramer-Rao bound at low 
SNR (i.e. cr^ is large), but suffers an estimation error floor at 
high SNR. Therefore, we go in for some improved estimators 
aiming at reducing this error floor 

III. Improved Estimators with Reduced Error 
Floor 

A. Method 1: Multiplying by Inverse of Matrix Q 

From the analysis in the last section, the estimation error 
floor comes from the fact that Q is not a perfect identity 
matrix. Therefore, a straightforward solution is to perform a 
linear transformation Q, such that fJQ = Iat. Since Q is known 
and always full rank for a given m-sequence, O = Q^. A 
new estimator is obtained by left multiplying the correlation- 
based estimator (|3]l by Q^: 



hi =Q-ih 



(9) 



which leads to an LS estimator that is in some extend similar to 
that proposed in [?]. The estimation error is = -i-Q^^Cw. 
The MSE of the estimator is: 
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As Q is a circulant matrix and all its elements are known for a 
given m-sequence, its inverse is easily obtained and does not 
need complex computations [?]: 

a b ■ ■ • b 

b a ■ ■ ■ b 



where 
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Replacing (fTTI ) into ( fTOb . the MSE becomes: 

2crl 



^hi N+1' 



(12) 



Comparing ( fT2] i and (|7]i, it can be found that the estimation 
error floor is removed by left multiplying by Q~^. In the 
meantime, the MSE is however approximately twice as much 
as the Cramer-Rao bound due to a noise power increase. 

B. Method 2: Multiplying by Inverse of Truncated Matrix Q 

Suppose the length of the CIR L is perfectly known, the 
estimate (O is truncated to L-length. The estimator becomes: 



h = Th = Qh 
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(13) 



where h is the L-length vector of the real CIR. T is an L x iV 
matrix whose left L x L submatrix is an identity matrix and 
rest parts are all O's, which represents deleting the last N — L 
elements of an A^-length vector. Q = TQT^ is an L x L 
circulant matrix that contains the first L rows and L columns 
of matrix Q and is still full rank. Similarly to (|9]l, we can left 
multiply estimator ( fT3] l by the inverse matrix Q to obtain: 

h2 =Q"^h = h+^Q"^TCw. (14) 

This new estimator is still unbiased. The related estimation 
error is = -^Q TCw which leads to an MSE given by: 

1^ /^r. ,H- 



= §Tr(Q-Q(Q-V).||Tr(Q-).(15) 

The inverse matrix Q ^ is the L x L circulant matrix with a 
similar form as (fTTI . The elements a and b are replaced by 



a = 1 + 



L- 1 



+ 2N - NL- L + 1- 



and 



N 



N'^ + 2N - NL- L + l' 
The MSE of the estimation is finally: 

N-L + 2 

^'^^ ^ N'^ + 2N -NL-L + l 



(16) 



Comparing (fTSl i and ( fT6] l, it can be found that the estimation 
MSE is reduced thanks to the truncation process. In the 
extreme case when N = L, the estimator ( fT4] i turns to (|9). 



When the channel length L is known, the Cramer-Rao 
bound of a training sequence based channel estimation is 
computed by replacing matrix P by P = PT"^ in ([8]): 



MS'£;>^Tr((p"p)-i~ 



{N-L + 2)al 



(17) 



m + 2N - NL- L + l 
This demonstrates that the proposed estimator (ITSt achieves 
the Cramer-Rao bound. 

C. Method 3: Subtracting Interference 

Another improved estimator can be obtained by subtracting 
the contribution of the interference from the correlation-based 
estimator. 

More precisely, suppose the CIR length L is known, we can 
rewrite the truncated CIR estimate (fTsl i as: 

1 
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(18) 



where A is the Lx L matrix that contains all the off-diagonal 
elements of Q and the main diagonal elements of A are O's, 
h is the L-length vector of the real CIR. Considering that 
A is known for a given m-sequence, we propose to use the 
estimated h to reduce the interference. An estimator is thus: 
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The estimation error is: = — Ah + -i-TCw — -^ATCw. 
RecaUing that Q = -^TCC^T^, the MSE of the estimate is: 

= iTr^AA(A)«+^(Q+AQA''-2K{QA''})) (21) 

where A is an L x L diagonal matrix with diagonal elements 
[af , ■ ■ ■ , aW, 3?{-} standing for the real part of a complex 
number. The first term of the MSE expression is: 

AA(A)^ = 
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The third term of the MSE expression is: 



^AQA^ 



Ar4 



The fourth term of the MSE expression is: 
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Fig. 2. MSB performance of different estimators with v = 420, N = 255, 
7V„„ = 165 in tlie TU-6 channel. 



TABLE I 

Complexities of the PN-based channel estimation methods. 



Estimator 


Complexity 


Circular con'elation-based estimator (|3] 


0(Af2) or 
O(Ar-logAf) 


Multiplying matrix inverse Af-length ^9) 
Multiplying matrix inverse L-length )14t 
Subtracting inteiference ||19I 


O(N^) 



Eventually, putting (|22li, (|23ll and (l24]i into (EB yields: 

_ N^ + {L-l){2-L-N) 2 {L-1){L^ -3L + 3) 

(25) 

Comparing (|25]l with (Q, in low SNR region, the MSE 
approaches the Cramer-Rao bound, while, in high SNR region, 
the estimation error floor is reduced by a ratio of approxi- 
mately (^)'. 

IV. Complexity Analysis 

In order to compare the computational complexity of each 
estimation method, we investigate the required basic oper- 
ations, i.e. multiplication, additions and FFT's etc. Table H] 
shows the complexity of the basic correlation-based estimator 
(O and the additional complexities required by the improved 
methods (|9|, (O and ( fT9] l. 

The A^-length circular convolution takes N multiplications 
and TV — 1 additions for each delay. That is N'^ multiplications 
and X (A^ — 1) additions for an A^-length CIR estimate. 
Therefore, the complexity of the circular correlation-based es- 
timator ^ is 0{N^). Considering that the circular convolution 
can also be computed by using the FFT, the computational 
complexity can be reduced to 0{N ■ log A^). 

As far as the improved estimators are concerned, additional 
complexities are needed by the estimation refinement process. 
For instance, in the estimator (|9]l, the matrix Q^^ depends only 
on the given PN sequence and can thus be pre-computed and 



Fig. 3. MSE performance of different estimators with u = 945, N = 511, 
N^„ = 434 in the TU-6 channel. 



TABLE II 

Profile of the TU-6 and HT channels 



Channel 


Tapl 


Tap2 


Tap3 


Tap4 


Taps 


Tap6 


TU 


Delay (/is) 





0.2 


0.5 


1.6 


2.3 


5.0 


Power (dB) 


-3 





-5 


-6 


-8 


-10 


HT 


Delay (/is) 





0.2 


0.4 


0.6 


15.0 


17.2 


Power (dB) 





-2 


-4 


-7 


-6 


-12 



stored 0. The matrix multiplication needs A^^ multiplications 
and A^ x (A^ — 1) additions. Therefore, the additional com- 
plexity is 0{N^). In the meantime, the matrix multiplication 
can be carried out in a reduced length L as done in estimator 
(fT4l) . For the same reason as in the pervious situation, the 
computation of Q does not need any additional effort. The 
corresponding additional complexity is therefore reduced to 
0{L^). For the estimator ([19}, the refinement of each channel 
tap needs {L — 1) additions. That is, L x (L — l) additions are 
required for the whole CIR estimate, which corresponds to a 
computational complexity of 0{L^). 

V. Simulation Results 

The simulation parameters are selected according to the 
specifications of the DTMB standard [?] where the sampling 
period is 1/7.56 /xs. The PN sequences are generated using 
the maximal linear feedback shift registers specified by the 
standard as well. The Typical Urban with six paths (TU-6) and 
Hilly Terrain (HT) channels specified in [?] are used in the 
evaluation. The power delay profiles of the two channels are 
given in Table [III The maximum delays of the TU-6 channel 
and HT channel are 5 /is and 17.2 /is which correspond to i = 
38 and L = 130 samples of the DTMB system, respectively. 
The performance of the different estimators is investigated 

'in fact, the matrix Q^' is quite structured. There are only two values - 
diagonal and off-diagonal elements. Hence, it is not necessary to store all the 
elements of the matrix. Instead, it is smarter to record the two values only. 
The cost of the storage is therefore negligible. 
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Fig. 4. MSB performance of different estimators with v = 420, N - 
Nf,„ = 165 in tlie HT channel. 



255, 



Fig. 5. MSE performance of different estimators with u - 
AT = 434 in the HT channel. 



945, N = 511, 



with different PN lengths and channels as shown from Figure |2] 
to Figure |5] Results are obtained from 1000 realizations of 
each channel. The correlation-based estimator corresponds to 
the classical approach found in the literature [?] while others 
are the improved versions introduced in section |III]A, B and 
C. From the figures, we can have the following observations: 

1 . Methods 'improved method 1 ' and 'improved metliod 2' com- 
pletely eliminate the estimation error floor, while method 
'improved metliod 3' reduces the error floor when the channel 
delay spread is shorter than the PN length. Multiplying the 
basic PN-correlation-based estimator by the inverse of the 
correlation matrix can create a perfect identity matrix which 
leads to estimators free of estimation error floor This can be 
seen from the performance of 'improved method 1 ' and 'improved 
method 2'. On the other hand, 'improved method 3' subtracts 
the interference components using the CIR estimates. It ap- 
proaches the Cramer-Rao bound in low SNR region, while the 
estimation error floor is reduced by a ratio of approximately 
{^ f in high SNR region. It indicates that the reduction is 
more notable when the channel delay is significantly small 
compared to the PN length. For instance, comparing Fig. |2] 
and Fig. IH 'improved method 3' obtains more improvement in 
the TU-6 channel than in the HT channel given the same PN. 
In addition, 'improved method 3' can achieve a lower estimation 
error floor when a longer PN is used (comparing Fig. Eland 
Fig. [51 or Fig. IH and Fig. [Sj. 

2. Method 'improved metliod V boosts the noise variance, 
while methods 'improved method 2' and 'improved method 3' do 
not. From this point, the latter two methods are the preferred 
ones because they outperform the classical correlation based 
approach whatever the SNR. In contrast, 'improved method 1' 
presents a performance back-off due to the noise component 
power boost, which leads to worse MSE results than the 
correlation based approach at low SNR. Last but not least, 
the estimator 'improved method 2' does not have estimation error 
floor and achieves the Cramer-Rao bound. 



VI. Conclusion 

In this paper, we have investigated the PN-correlation- 
based channel estimator for TDS-OFDM. Aiming at reducing 
the estimation error floor encountered by the classical PN- 
correlation-based estimator, we have proposed three improved 
estimators which exploit the correlation property of the m- 
sequence and the knowledge of the channel length. It has been 
shown through mathematical derivations and simulations that 
the new proposed estimators approaches or even achieve the 
Cramer-Rao bound. 

Acknowledgment 

The authors would like to thank the European CELTIC 
project "ENGINES" for its support of this work. 



